High cumulants from the 3-dimensional 0(1,2,4) spin models 
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Considering different universality classes of the QCD phase transitions, we perform the Monte 
Carlo simulations of the 3-dimensional 0(f,2,4) models at vanishing and non-vanishing external 
field, respectively. Interesting high cumulants of the order parameter and energy from O(l) (Ising) 
spin model, and the cumulants of the energy from 0(2) and 0(4) spin models are presented. The 
generic features of the cumulants are discussed. They are instructive to the high cumulants of the 
net baryon number in the QCD phase transitions. 
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I. INTRODUCTION 

A thermodynamic system may undergo a phase tran- 
sition when temperature and/or pressure vary. For the 
strong interacting QCD system, it is known that the 
hadronic matter will change to partonic phase at very 
high temperature and density [l|. The data from ultra- 
relativistic nuclear-nuclear collisions at RHIC has shown 
that the partonic phase of QCD matter-quark-gluon 
plasma (QGP) has been formed 0]. One of the main 
goals of current ultra-relativistic heavy-ion collision ex- 
periments is to map the QCD phase diagram Q. 

The critical point at the QCD phase boundary is par- 
ticular interested, as the large fluctuations are expected. 
The sensitive probes of QCD critical point are the high 
cumulants of the conserved charges, i.e., net baryon num- 
ber, net electric charge, and net strangeness (j-|9(. 

The net baryon number fluctuations have been studied 
using lattice QCD simulations and QCD effective mod- 
els fl0l - tl3j . However, due to the difficulties of the lattice 
calculations and model estimations, the hi gh c umulants 
of the net baryon number are still not final (Til. HH ] . The 
universality of the critical fluctuations allows us to study 
the relevant cumulants in a simple system, where the high 
cumulants can be precisely obtained. 

It is known that the QCD critical point, which termi- 
nates the first order phase transition line, is in the same 
universality class as the 3-dimensional Ising model [iH - 
[j~9j . So the high cumulants of order parameter and energy 
in the vicinity of the critical point of the 3-dimcnsional 
Ising model arc called for. They should be a good refer- 
ence for various relevant calculations. 

In the chiral limit, the chiral phase transition for 2- 
flavor QCD belongs to the same universality class as the 
3-dimensional 0(4) spin model The singular behav- 
ior of the net baryon number fluctuations is expected to 
be governed by the universal 0(4) symmetry group [2(ij |. 
Owing to the staggered fermions in lattice calculations, 
the (2+l)-flavor chiral phase transition may fall into the 
3-dimcnsional 0(2) universality class [2lT - l23j . 

In this paper, we present the high cumulants of the 
order parameter and energy in the vicinity of the critical 



temperature from the 3-dimcnsional Ising model, and en- 
ergy from the 3-dimcnsional 0(2) and 0(4) spin models, 
respectively. The simulation is performed in a finite-size 
system. The obtained results should be instructive to a 
finite system formed in ultra-relativistic heavy-ion colli- 
sions. 

The paper is organized as follows, firstly, the cumu- 
lants of the order parameter and energy from the O(N) 
spin models are derived in section II. Then, their rela- 
tions to the net baryon number fluctuations in the QCD 
phase transitions are discussed. In section III, the high 
cumulants from the 3-dimensional Ising, 0(2) and 0(4) 
spin models are presented and discussed. Finally, the 
generic features of the cumulants from the three models 
are summarized in section IV. 



II. CUMULANTS IN THE O(N) SPIN MODELS 

The 0(iV)-invariant nonlinear er-models (O(N) spin 
models) are defined as, 



pn = - J Y & ■ S - H ■ Y Si 



(1) 



where T-L is the Hamiltonian, J is an interaction energy 
between nearest-neighbor spins and H is the ex- 

ternal magnetic field. J and H are both reduced quan- 
tities which already contain a factor ft = 1/T. Si is a 
unit vector of iV-components at site i of a d-dimensional 
hyper-cubic lattice. It is usually decomposed into the 
longitudinal (parallel to the magnetic field H) and the 
transverse component 



Sf e H + S, 



(2) 



where e# — H/H. For the 3-dimensional Ising, 0(2), 
and 0(4) spin models, d = 3, and N = 1, 2 and 4, re- 
spectively. 

The (reduced) free energy per unit volume is 



f(T,H) = ~]nZ, 



(3) 
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where Z = J]J l d N S t S(Sf - 1) exp(-/3^ + HVS^ ) is the 
partition function. E = — j\ Si ■ Sj is the energy of a 

spin configuration, 5 1 " = y J2 i S 1 ]' is the lattice average of 
the longitudinal spin components, V = L 3 is the volume, 
and L is the number of lattice points of each direction. 

As we known, the cumulants of the order parameter 
are the derivatives of the free energy density (Eq. (O) 
with respect to H. We can get the cumulants of the 
order parameter from the generating function |24j |. 



ln(e 
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(4) 



x=0 



So the first, second, third, fourth and sixth order cumu- 
lants of the order parameter arc as follows, 



4 = (55ll 4 )-3(^H 2 ) 2 , 
4 = (^H 6 )-10{^ll 3 ) 2 + 30( ( 55ll 2 ) 3 
-15<<5S I|4 )(5S 112 ), 



(8SU ) 



(5) 



where tfS" = S'l - (S 11 ), and ftf is the magnetization 
(order parameter) of the system. At vanishing external 
magnetic field, due to the spatial rotation symmetry of 
the O(N) groups, such defined order parameter is zero. 
In the case, an approximated order parameter definition 
is suggested as, M = (|i Ei-^l) H- 

On the other hand, the cumulants of the energy are 
the derivatives of the free energy density with respect to 
the temperature T. The generating function is 



3— \n(e xE ) 
dx 11 x ' 



(6) 



x=0 



So the first, second, third, fourth and sixth order cumu- 
lants of the energy are as follows, 



--(E), ^ = (5E 2 ), 
-- (SE 4 ) - 3(SE 2 ) 2 , 
-- (5E 6 ) - 10{5E 3 ) 2 + 
- 15(5E 4 )(5E 2 ), 



Ki = (SE 3 }, 
30(5E 2 ) 3 



(7) 



where SE = E - (E). 

In the vicinity of the critical point, the free energy 
density (Eq.©) can be decomposed into two parts, the 
regular and singular parts. The critical related fluctu- 
ations are determined by the singular part. It has the 
scaling form 



f s (t,h) = r d f s (iy%iy*h). 



(8) 



Here t = (T — T c )/T and h = H/H a are reduced tem- 
perature and magnetic field, To and Hq are the normal- 
ized parameters. T c is the critical temperature. y t and 
Uh are universal critical exponents. In our simulation, 



we set J = (3 and choose the approximate critical tem- 
peratures, T c Isin s = 4.51 [Hi, T c ° (2) = 2.202 [H, and 
T c ° (4) = 1.068 [11 for the 3-dimensional Ising, 0(2) and 
0(4) spin models, respectively. 

In order to map the result of the 3-dimcnsional Ising 
model to that of the QCD, the following linear ansatz is 
suggested 

t^T-T cp + a(n- n cp ), h « fi - fi cp + b(T - T cp ). (9) 

T cp and fi cp are the temperature and chemical poten- 
tial at the QCD critical point, respectively, a and b are 
two parameters determined by QCD. The baryon-baryon 
correlation length diverges with the exponent y t and ex- 
ponent yh when the critical point is approached along the 
i-direction and ^.-direction, respectively flit ] . 

The cumulants of the net baryon number are the 
derivatives of the QCD free energy density with respect 
to fi. In the vicinity of the critical point, it is the com- 
bination of the derivatives with respect to t and h in the 
3-dimensional Ising model. Since yh is larger than y t [30| . 
the critical behavior of the net baryon number fluctua- 
tions is mainly controlled by the derivatives with respect 
to h, i.e., the fluctuations of the order parameter in the 
3-dimcnsional Ising model. 

The singular part of the free energy density for the 
chiral phase transition is suggested as Q 



fs(T,fl q ,h) 

2^4 



AhW>ff(z), z = t/h? s , 



(10) 



where (3 and 6 are the universal critical exponents of the 
3-dimensional 0(4) spin model. ff(z) is the scaling func- 
tion. The reduced temperature t and external field h are 
expressed as follows 



t 



1 m q 



111) 



Here T c is the critical temperature in the chiral limit, 
is a parameter determined by QCD [23| • The net baryon 
number susceptibility is the derivative of free energy den- 
sity with respective to the chemical potential fi q . From 
Eqs. (jTU)) and (|TT1) . wc can get the derivatives of the free 
energy density with respect to T, to chemical potential 
fiq = fJ,q/T at £i q = 0, and to ft q ^ respectively, 



df/T 4 



A 



dT r, 



{t T c y 



-hP- a -*V fls f} l \z), (12) 



df/T 4 



0£? 
df/T 4 



-A(2K q )^ 2 h^ 2 - a - n ^^ s & /2 \z), 



fj. q =0 



where n is even in Eq. ([1 



(13) 



- A (2 Kq n^rh( 2 - a - n y^f^(z), 



(14) 
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Comparing Eqs. ([13]) and {H) with Eq. ((H), we can 
see that the scaling form of derivatives of the free energy 
density with respect to T and [i q are equivalent. So the 
net baryon number fluctuations are related to the deriva- 
tives of the free energy density with respect to T in the 
3-dimcnsional 0(4) spin model. The n-th order cumulant 
of the energy from the 3-dimcnsional 0(4) spin model is 
relevant to the 2n-th (or n-th ) order cumulant of the net 
baryon number at fj. q = (or [i q ^ ) in the chiral phase 
transition. 

The cumulants of the order parameter from the 3- 
dimcnsional 0(2) and 0(4) spin models are the deriva- 
tives of the free energy density with respect to the exter- 
nal field. Their critical behavior has been presented and 
discussed in Ref. 1311 . 



III. CRITICAL BEHAVIOR OF THE HIGH 
CUMULANTS 

The Monte Carlo simulations of the 3-dimensional 
Ising, 0(2), and 0(4) spin models are performed by the 
Wolff algorithm with helical boundary conditions [32T |. 
We start simulation without magnetic field, and then 
modify the algorithm to include a magnetic field H = 
0.05 and fl = 0.1@. 

As we know, the valid region of system size of the finite- 
size scaling varies with magnetic field and observable. 
The typical size of an observable at a given magnetic 
field is determined by the saturation of size dependence, 
as shown in Ref. 



34 



The typical system sizes for each 
kind of cumulants at a given magnetic field and model 
arc listed in Table I. 



TABLE I: The typical system size for each kind of cumulant 
at a given external field and model 



■----^observable 
field(H)~~~~----^_ 


4(0(1)) 


«f(o(i)) 


ref(0(2)) 


ref(0(4)) 





24 


20 


20 


20 


0.05 


12 


10 


10 


8 


0.1 


8 


8 


8 


8 



In order to compare the basic structure of the cumu- 
lants with and without external fields, each cumulant of 
the energy at vanishing and non- vanishing external fields 
is plotted in an identical sub-figure, and rescalcd to unity 
by its maximum or minimum (except for the first order 
cumulant of the energy from 0(2) and 0(4) spin models), 
as shown in Fig. 1, Fig. 4, and Fig. 5. For the Ising model, 
the cumulants of the order parameter at vanishing and 
non- vanishing external fields arc quite different, which 
may be caused by the definitions as discussed in section 
II below Eq. ([5]). So their cumulants at non- vanishing 
and vanishing external fields arc presented in Fig. 2 and 
Fig. 3, respectively. 

The cumulants of the energy, i.e., ref, re^f, ref and , 
from the 3-dimensional Ising model at H = 0.1,0.05,0 
are shown in Fig. 1(a) to 1(d), respectively. We can see 



that the basic features of the cumulants, i.e., the patterns 
of the fluctuations, are not influenced by the magnitude 
of the external field. With the appearing or increasing of 
the external field, the whole critical region is amplified 
and shifted toward the high temperature side. 

In the vicinity of the critical temperature, ref has 
a maximum peak, ref" oscillates and its sign changes 
from positive to negative when temperature increases and 
passes the critical one. ref has two positive maximums lo- 
cating at the two sides of T c , respectively. The minimum 
between them is negative. In contrast to the two pos- 
itive maximums of ref, ref" has two negative minimums 
and one positive maximum at the critical temperature 
region. The sign change for the cumulants of the energy 
starts at the third cumulant. It happens twice in the 
fourth and sixth order cumulants. 

The cumulants of the order parameter at non- vanishing 
external fields, i.e., ref, ref, ref and ref, from the 3- 
dimcnsional Ising model are presented in Fig. 2(a) to 
2(d), respectively. We can see that the influences of the 
external field are similar to those as discussed in Fig. 1. 

In the vicinity of T c , ref shows the same peak structure 
as that for the energy, ref has a negative valley and no 
sign change in the critical region, ref shows a obvious pos- 
itive maximum and a very small negative minimum when 
the temperature increases and passes the critical one. ref 
oscillates from positive to negative, and the negative val- 
ley is more obvious than that in ref. Here the sign change 
starts at the fourth order cumulant in Fig. 2(c). 

The high cumulants of the order parameter from the 
3-dimcnsional Ising model at vanishing external field are 
shown in Fig. 3(a) to 3(b), respectively, ref in Fig. 3(a) 
shows a narrow and sharp peak at the critical temper- 
ature region. From Fig. 3(b), we can see that, both ref 
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FIG. 1: (Color online), ref (a), ref (b), ref (c) and ref (d) at 
H — 0.1, 0.05, from the 3-dimensional Ising model. 
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and nf oscillate, but the former changes from negative to 
positive, while the latter changes from positive to nega- 
tive with the increasing temperature. The generic struc- 
ture of Kg is similar to that of nf, having two positive 
maximum locating at the two sides of T c and a negative 
minimum between them. So the generic structure of the 
high cumulants at non-vanishing external field is quite 
different from that at vanishing external field. The sign 
change in the former case appears in the fourth order 
cumulant, while the third one in the later case. 

Comparing the energy fluctuations in Fig. 1 with the 
order parameter fluctuations in Fig. 2 and Fig. 3, we can 
see that the generic structure of the same order cumulant 
of the energy is different from that of the order parame- 
ter, except for the second order one. The way of the sign 
change of n§ at H = as shown in Fig. 3(b) is consistent 
with the expectation from the effective model [3]. The 
generic structure of ref at non-vanishing external field 
as shown in Fig. 2(c) is very similar to that from the 
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FIG. 3: (Color online). «f (a), «f , k| and k§ (b) at H = 
from the 3-dimensional Ising model. 
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The cumulants of the energy, i.e., ttf , /cf 7 , k§ and 
nf , from the 3-dimensional 0(4) spin model at H = 
0.1, 0.05, are presented in Fig. 4(a) to 4(d), respectively. 
Again, the external field shows the similar influences as 
discussed above, nf increases with the temperature. K% 
has a peak and is positive in the whole critical temper- 
ature region, ref" oscillates and changes from positive to 
negative with the increasing temperature, nf has two 
maximums and a minimum between them. The behavior 
of k,2 > K 3 an d ftf is similar to that of the energy fluctu- 
ations from the 3-dimensional Ising model, cf. Fig. 1(a) 
to 1(c). 

As we discussed in section II, at vanishing baryon 
chemical potential, xi > Xi > xi an d xi m the chiral 
phase transition are corresponding to nf , , and 
nf from the 3-dimcnsional 0(4) spin model. The positive 
peak of k§ is consistent with xf from the calculations 
of lattice QCD and the estimations of the PNJL [lj] 
and PQM models The sign change of is also 

observed in xi from the PQM model [aHl- 

The cumulants of the energy, i.e., tzf, nf, «f and 
nf , from the 3-dimensional 0(2) spin model at H = 
0.1, 0.05, are presented in Fig. 5(a) to 5(d), respectively. 
We can see that each sub-figure in Fig. 5 is similar to that 
in Fig. 4. Each order cumulant of the energy from the 
3-dimcnsional 0(2) and 0(4) spin models shows qualita- 
tively similar temperature dependence in the vicinity of 
the critical temperature. 
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FIG. 5: (Color online), nf (a), «f (b), k§ (c) and >tf (d) at 
H = 0.1, 0.05, from the 3-dimensional 0(2) spin model. 



IV. SUMMARY 

In this paper we perform the simulations of the 3- 
dimensional Ising, 0(2) and 0(4) spin models at a 
given system size at three different external fields H = 
0.1,0.05,0. The critical behavior of the high cumulants 
of the order parameter and energy in the 3-dimcnsional 
Ising model, and the cumulants of the energy in the 3- 
dimensional 0(2) and 0(4) spin models is presented, re- 
spectively. We find that the external field does not influ- 
ence the generic structure of the cumulants, except the 
cumulants of the order parameter from the 3-dimensional 
Ising model. But it widens the temperature region of the 



critical fluctuations. 

For the 3-dimcnsional Ising model, the generic struc- 
ture of the high cumulants of the energy are different from 
that of the order parameter. For the energy fluctuations, 
the sign change starts at the third order cumulant no 
matter with or without external field. So docs the order 
parameter at vanishing external field. At non-vanishing 
external field, the first sign change of the order param- 
eter fluctuations appears at the fourth order cumulant. 
The common feature is that the higher the order of the 
cumulant, the more complicated the fluctuation pattern 
is. 

The critical behavior of the third order cumulant of 
the order parameter at vanishing external field, and the 
fourth order cumulant of the order parameter arc consis- 
tent with the expectations of the effective model and the 
CT-model, respectively. 

For the 3-dimensional 0(2) and 0(4) spin models, the 
behavior of the second to fourth order cumulants of the 
energy is similar to that from the 3-dimensional Ising 
model. The sign change also starts at the third order 
cumulant. 

The net baryon number fluctuations based on the 0(4) 
spin model are qualitatively consistent with the calcula- 
tions from the lattice QCD, and expectations from the 
PNJL and the PQM models. Our results also show that 
at vanishing chemical potential, the sixth order cumulant 
of the net baryon number is necessary in order to observe 
a sign change in the chiral phase transition. 
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